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I. INTRODUCTION 



A quantum operator (or quantum superoperator) $ on the collection Ai^ of complex N x 
matrices is a completely positive, trace preserving linear map. The quantum operator $ is 
unital provided that $(/) = /, that is, provided that $ fixes the identity matrix. A density 
matrix, which represents the state of a quantum system, is a positive matrix (Hermitian 
with nonnegative eigenvalues) having trace one. The properties of density matrices are 
thus preserved via the action of a quantum operation. Of course, density matrices are 
mapped to density matrices under any trace-preserving positive superoperator. (A positive 
superoperator, by definition, takes positive matrices to positive matrices.) The requirement 
that a quantum operation be completely positive rather than simply positive is based on the 
viewpoint that $ represents the "restriction" of a positive operator on a larger system. 
By definition, a superoperator $ on A^tv is completely positive provided that / ® $ : 
M-m ® M.N ^ ® A^Af is positive for all positive integers m. 

A Density matrix p G M.2 represents the state of a two-level quantum system — a one 
qubit system. It's not difficult to show that such matrices have the following "Bloch" 
representation: 

P = ^i%^. (1) 



where (ri,r2,r3) belongs to the closed unit ball of M'^ and where 
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,0-3 = 




1 




i 




-1 



0-1 



are the usual Pauli matrices. The vector r = {ri,r2,rs) appearing in is called the Bloch 
vector or coherence vector of p. The correspondence between elements of the closed unit 
ball of and density matrices is complete for two-level systems; that is, a trace-one matrix 
p E M.2 is positive if and only if it has representation where |r| < 1. 
A density matrix p E M.n has a representation analogous to (fT]): 

Af2-1 



p 



1 + 



N{N 



^ E r.A. 



i=l 



where now the Bloch vector r* belongs to the closed unit ball of . 



(2) 



and {Xi}i=i consists 



of elements of at having the following properties: 

Xi is self-adjoint (AJ = Aj), trAj = 0, and tr(AjAj) 
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26. 



for two-level systems, the collection of Bloch vectors r from (0) that corresponc 
matrices is a proper subset of the unit ball of M^^"^, recently characterized in 0, 



One may, for example, take A to consist of generators of SU(A^) (see j2; 1^ In our 
representation of density matrices, we have adopted the normalization factor {N{N — 
l)/2)^/^ found in P|, which forces a pure-state density matrix to have a Bloch vector r of 
norm 1. Note that together with the identity matrix / constitutes an orthogonal 

basis oi M.]^ relative to the Hilbert-Schmidt inner product. In contrast to the situation 

to density 
5|. We will 

be concerned with n-qubit systems, which means that in 0, = 2" and each element of 
{Aj}^]"^ may be taken to be an appropriately normalized tensor product of n Pauli matrices, 
where ctq := / is included in the "Pauli collection" but not all factors in the product can be 
ctq — see Section El below for details. 

Because unital quantum operators are completely positive and preserve both the trace 
and the identity, associated with any such operator $ on A^jv? there is an (A^^ — 1) x (A^^ — 1) 
real matrix M$ such that 



where c = N(N — l)/2 is the normalizing constant. We call M$ the Bloch matrix of $. 
For N = 2, Bloch matrices M that correspond to unital quantum operators are charac- 



terized in 



ll (see also m and 



a 



12| |) in terms of signed singular values of M: 



^ — ^ \di — and 1 + c?3 > |c/i + ^2!, (3) 

where di, d2, and d^ may be taken to be the singular values of M if det M > and may be 

det M < 0. (For information 
10lli3l, e.g.) The inequalities Q 



taken to be the additive inverses of the singular values of M i 



about singular values of matrices, the reader may consult 
are the Algoet-Fujiwara conditions for a Bloch matrix corresponding to a unital quantum 
operator on the Bloch ball jl| . In Section 4 of this paper, we present a short and completely 
elementary derivation of the Algoet-Fujiwara conditions, and in Section 5 we show how our 
methods yield a description of "diagonal" quantum operations on n-qubit systems, obtaining 

ihe n-qubit analogue of the Algoet-Fujiwara conditions (see Theorem IV. II below) . C. King 
has obtained a related description of diagonal quantum operators on three-state systems 

i.e., single qutrit systems). 
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In the next section, we discuss operator-sum representations of superoperators that map 
Hermitian matrices to Hermitian matrices. As is well-known any superoperator $ on 
Aiisr that preserves Hermiticity must have the form 

k 

i=i 

where for each j, Aj is a.n N x N matrix and ej G {1,-1}. We show in Proposition III. II 
below that if $ is completely positive and the operator elements {Aj} for $ are linearly 
independent, then ej = 1 for j = 1,2, ... ,k. Proposition 111.11 is the key lemma in our work 
of Sections 4 and 5. 

Before concluding this introduction, we should add a remark about non-unital quantum 
operators. These superoperators correspond to affine mappings: {I + cf ■ X)/N {I + 
c{Mf + t) ■ A) /N. For = 2^ffine mappings r — > Mr + t corresponding to quantum 
operators are characterized in (l2l |. 



II. SIGN PATTERNS IN OPERATOR-SUM DECOMPOSITIONS 

Let $ be a superoperator ow M.^ and suppose that for some positive integer k there exist 
N X N matrices Ai, A2, . . . , Ak along with "signs" Ej G {—1,1} such that 

k 

Hp) = J2^,A,pA]. (4) 

i=i 

The expression of the right of (@)) is called an operator-sum decomposition of $ and {Aj}'^^^ 
is corresponding set of operator elements. Operator-sum decompositions in which Sj = 1 

n nn 

for every j model system-environment interactions ([l^; see also la nil)- For this reason, 
operator elements are sometimes called "interaction operators" . 

Observe that if $ has an operator-sum decomposition (ji}, then $ preserves Hermiticity; 
that is, $(p)^ = $(p) whenever p is Hermitian. In de Pillis shows that every superop- 
erator on A^Tv that preserves Hermiticity has an operator-sum decomposition. For example, 
by de Pillis's result, the transpose operator $t on M.2 defined by $t(p) = must have 
an operator-sum decomposition. A simple calculation shows that one such decomposition is 
given by 

^ r \ ^ I (Ji ai Cr2 0-2 (73 (73 , . 



Operator-sum representations are highly non-unique; for instance $t is also given by 



$t(p) = 
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It's obvious that $r is a positive superoperator, preserving both Hermiticity and eigenvalues. 
On the other hand $t is the canonical example of a positive operator that is not completely 
positive. Can the fact that is not completely positive be deduced from the presence of the 
negative sign in the operator-sum decompositions for displayed above? Proposition lII.il 
below shows that the answer to this question is yes. This is not an entirely trivial matter. 
For example, the identity superoperator operator $/, which is obviously completely positive, 
is given by 

$,(P) = (V2/)p(V2/) - Ipl. (6) 

As we discussed in the Introduction, quantum operations are completely positive. Ob- 
serve that since Hermitian matrices are differences of positive matrices (immediate from the 
spectral decomposition), positive (and hence completely positive) superoperators must pre- 
serve Hermiticity. Thus any quantum operator $ has an operator-sum decomposition (^. 
In Theorem 1 of jsl , Choi shows that a completely positive operator $ has an operator-sum 
decomposition Q in which each sign is positive {Sj = 1 for every j). Of course this doesn't 
mean that every operator-sum decomposition of a completely positive map must feature 
only positive signs, as © shows. Choi does not prove his Theorem 1 as a corollary of de 
Pillis's theorem for Hermiticity-preserving superoperators. Rather, he gives an elegant in- 

n 

dependent proof that also yields de Pillis's characterization (as Choi points out 3, p. 277]). 
More important for our purposes is that Choi investigates the relationship between differ- 
ent operator-sum representations of the same superoperator, proving 0, Remark 4] that if 
{Aj}^^^ and {Ej}^^^ are collections of operator elements for the same superoperator on M.^ 
and if {Aj}^^^ is linearly independent in /iAn, then there is an isometric mx k matrix [ajn] 
such that for each j G {1, 2, . . . , m} 

k 



(7) 



n=l 



Note that because [aj-n\ is an isometry, we must have m > fc; if both {Aj}^^^ and {Ej}'^^^ 
are linearly independent collections of operator elements for the same superoperator, then 
Choi shows m = k and the matrix of scalars [ctj^] relating them by ((7j) is unitary. 



Proposition II. 1 Suppose that ^ : M.^ ^ A4n is given by 

k 

i=i 

where {Aj}'^^^ is linearly independent in M.j~q and Sj G {—1, 1} for j = 1,2, . . . , k. Then 
$ is completely positive if and only if Sj = 1 for each j G {1,2, . . . , k} . 

Proof. If each sign "ej" is positive, then $ is easily seen to be completely positive (and the 
independence of {Aj} is irrelevant). 

Conversely, suppose that $ is completely positive and has the form displayed in the 
statement of the proposition with {Aj}^^^ independent. We assume, in order to obtain a 
contradiction, that some of the signs ej are —1. Without loss of generality, we take Sj = —1 
for j = 1 to p for some p E {1,2, . . . ,k — 1}. (Clearly, not all of the signs can be negative: 
the linear independence of the set of operator elements means that no element Aj is the zero 
matrix so that if all signs were negative, $ would map positive matrices to negative ones 
and hence couldn't be completely positive). 

Because $ is completely positive, Choi's work shows that there exists an operator-sum 
decomposition for $ with all signs positive: 

m 

i=i 

We have for every N x N matrix p, 

k m 

or 

k m p 

E ^^/^^J = E ^^p^^ + E "^M^ (8) 
j=p+i j=i j=i 

so that we have two different operator-sum representations for the same superoperator p i-^ 
Yl'j=p+i AjP^]- Thus, in particular, there are scalars (q;„) (forming one row of the isometric 
matrix relating the operator elements on the left of (jHj) to those on the right) such that 

k 

^1 = E 

n=p+l 

contradicting the linear independence of {Aj}*'^^. QED 

In the sequel, we will use the following quick corollary of the preceding proposition. 
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Corollary II. 2 Suppose that ^ is a completely positive superoperator on A4n having the 
representation 

k 

i=i 

where {Aj}^-^^ is linearly independent in J^jy and [3j is real for j = 1, 2, . . . , k. Then 
(3,>0forje{l,2,...,k}. 

Because the Pauli matrices ai, a2, and together with the 2x2 identity matrix form 
a linearly independent set in 7Vl2; the preceding corollary shows that the superoperator $ 
defined on by 

$(p) = Poaopao + Piaipai + p2(^2pcr2 + Pscrspc^s (9) 

is completely positive only when j3j > for j = 0, 1, 2, 3. It's not difficult to obtain a 
characterization of positivity for the superoperator $ defined by Q. The characterization, 
presented in the next proposition, shows that if $ is positive but not completely positive 
and $ is written in the form of above, then exactly one of the scalars f3j will be negative, 
as illustrated in equation 

Proposition II. 3 The superoperator $ : -^2 defined by (0j is positive if and only 

if every pair from {(3j}^^Q sums to a nonnegative number: 

Po + Pi> 0,Po + /32 > 0,/3o + /53 > 0,/5i +P2>0,Pi+P3> 0,/32 + /^a > 0. (10) 

Proof. The following observation will facilitate some calculations in the proof; it will also 
play in a crucial role in the final two sections of this paper. 

(7iajai = ±aj, i, j G {0, 1, 2, 3}, (11) 

where the sign is positive when z = 0, j = 0, or i = j, and negative otherwise. 

We assume that $, defined by is positive and obtain the inequalities stated in the 
proposition. Suppose p is a positive matrix. Without loss of generality we will assume that 
p has trace 1 and hence has the form (/ + Viai + r20"2 + r3cr3)/2, where r = (ri,r2,r3) lies 
in the unit ball of M'^. A simple calculation shows 

Sol + nsiai + r2S2(T2 + r^ssas . 
Hp) = ^ , (12) 



where sq = /3o + A + /32 + /^s, si = + f^i - f^2 - Ps, S2 = (3o - (3i + (^2 - Ps, and 
S3 = Po — Pi — P2 + Ps Because $(p) is positive, its trace is nonnegative; thus, 



tr($(p)) = So > 0. 

If So = 0, then the positive matrix is the zero matrix (independent of p), which, in 
view of ()12j). forces Si, S2, and S3 to be zero as well. It follows that Pj = for each j, and 
the inequahties (|Tn|l hold for this case. 

Suppose tr($(p)) > 0. Then we may rewrite the right-hand side of (|T^ as follows: 

So " ^ " ^ j , (13) 

which represents a positive matrix if and only if 

Sl\ , ( S2\ , f S3 



Substituting, respectively, r = (1,0,0), (0,1,0), and (0,0,1) into (fT^ yields the following 
three inequalities 

< So, |S2| < So, |S3| < So, 

which, in turn, yield the desired inequalities (fTUI) . 

For the proof of the converse, suppose that $ acts on the positive matrix (J + r-a)/2) and 
that the inequalities (fTIH) hold. Adding the first and last inequalities of (jlOj) . we must have 
Po + Pi + P2 + P3 ^ 0- Po + Pi + P2 + P3 = 0, then by grouping summands appropriately, 
one obtains that each of the inequalities of (fTUI) must be an equality and it follows from ()12|) 
that $ is the zero operator. On the other hand, if Po + Pi + P2 + Ps > , then it's easy to 
see that the quotients multiplying ri, r2, and in (fT^ must each have absolute value less 
than or equal to 1 and this shows that the quantity on the left of (jl4|) is < Ti + r| + r|, 
which is < 1 since p is positive. Thus $(p) is positive, as desired. QED 



III. UNITARY SUPEROPERATORS AND ROTATIONS OF THE BLOCK RE- 
GION 

In this section, we summarize known information about quantum operators having an 
operator-sum decomposition with single unitary operator element. 



8 



Suppose that $ : M.n M.n is unitary in the sense that it has an operator-sum 
representation of the form 

$(p) = UpU\ 

where ?7 is a unitary N x N matrix. Clearly such unitary $'s are unital quantum operators. 

It's well known (see, e.g., Exercise 8.13]) that if $ is unitary and acts on A^2 (the one 
qubit situation), then its Bloch matrix M<i, is a rotation matrix on M.^; that is, M^M,^ = I 
and det(M$) = 1. Furthermore, it's not difficult to show that the correspondence between 
rotation matrices on and unitary superoperators is complete in the N = 2 setting; that 
is, given any rotation matrix A there is a unitary 2x2 matrix U such that 



2 / V 2 

Now suppose N > 2 and ^ : A4n ^ Aijy is a unitary quantum operator. As one would 
expect, once again $ has a Bloch matrix M$ that is a "rotation", where by rotation we 
mean M$ is orthogonal (MjM$ = /) and orientation preserving (det M = 1). It's very easy 
to see that M$ must be orthogonal: let ^(p) = UpW and p = (/ + cr ■ A)/A^ and note 

^ + (1 - ^)|rT = tr(p^) = tr(f/pf/tf/pf/t) = tr ^L±^^L:A^ = ^ + " 

(15) 

so that |M$r| = |r|. Thus is an isometry and since it has real entries, M$ is orthogonal. 
The proof that det(M$) = 1, which we will also present, requires a bit more effort. Because 
U is unitary, there is an orthonormal basis {\vj))jLi of consisting of eigenvectors of U 
with corresponding eigenvalues (e*'*^ )^^, where the Cj's are real. For s G [0, 1], define 

j 

so that Uq = U and Ui = I. Let : Mn ^ Mn he given by $s(p) = UspUl and let 
be the corresponding Bloch matrix. We have already shown that is orthogonal for each 
s in [0, 1]. Hence det(M$J = ±1 for each such s. Since det(M$J varies continuously with 
s and since det(M$J = det / = 1, we see det(M,j,) = det{M^^) = 1, as desired. 

When N > 2, the correspondence between rotations and unitary quantum operators is 
complicated; for example, the angle 6 between pure-state Bloch vectors fl and must 
satisfy cos(6') > — 1/(A^ — 1), or, equivalently, rl ■ > — 1/(A^ — 1) 



Returning to the one-qubit situation, suppose that M is an arbitrary 3x3 real matrix 
and the hnear superoperator $m : M.2 M2 is defined by 

^ fl + f-a\ I + Mr-a 

^Mi^^]= ^ . (16) 



An interesting problem is to determine when $m is a quantum operation. An obvious 
necessary condition is ||M|| < 1, where ||M|| = max{|M?^ : |?^ = 1}. A complete descri ptio n 
of those M such that is a quantum operator may be found in |l] (see also and 
This description is based on the singular-value decomposition of M, which, in turn, easily 
yields the following. 



Proposition III.l Suppose M is a 3 x 3 matrix with real entries. Then there exist 3x3 
rotation matrices A and B as well as a 3 x 3 diagonal matrix D with real entries such that 

M = BDA. 

Moreover, if det M > 0, then the diagonal entries of D may he chosen to he the singular 
values of M listed in decreasing order; otherwise, the diagonal entries of D may he chosen 
to he the negatives of the singular values of M listed in increasing order. 

Let $jvf be the unital superoperator (|T6|l on M.2 induced by the real 3x3 matrix M. Let 



M = BDA be the factorization of M promised by Proposition IIII. 1| and let let Ua and Ub 
be the unitary 2x2 matrices such that 

Finally, let be the unital superoperator defined by 

, I + f ■ a\ I + Dr ■ a 

Note that 



2 / 2 



<I>m(p) = (^o$^o1])(p), (18) 

where \I/(p) = UbpU'^b ^^'^ ^(p) ~ UapU\. Because \1/ and Vt and their inverses are quantum 
operations and compositions of quantum operations are quantum operations (jl8j) shows that 
is a quantum operation if and only if is a quantum operation. Thus, to characterize 
unital quantum operators on the Bloch ball, one need only understand which diagonal 
matrices M are such that $m is a quantum operator. Necessary and sufficient conditions on 
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the diagonal entries of M (which are given by Q in the Introduction) that ensure M induces 
a_ quantum operation are obtained in jl| as well as 1^ and The method employed in 



Q and 



to obtain the conditions is based on the proof of Theorem 



1 of 0: 



one analyzes 



the positivity of 



E 







where / is the identity on A^2 and where E is the 4x4 matrix composed 4 elementary 2x2 
blocks: 

■"l 1 
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The method employed in [7[ involves starting with an operator-sum decomposition of the 
unital quantum operator in question and expressing its operator elements as linear combi- 
nations of the Pauli matrices. 

In the next section, we take a new approach to characterizing the diagonal matrices 
corresponding to quantum operators on the Bloch ball. Our approach is based on our 
work with sign patterns in operator-sum decompositions in Section |n] and allows convenient 
generalization to the n-qubit situation. 



IV. DIAGONAL QUANTUM OPERATORS ON THE BLOCH BALL 



Suppose that 

dii 

D= ^22 
^33 

and $D is the linear superoperator on A^2 defined by 

I + r ■ (y\ I + Df ■ a 



D 



2 / 2 (1^) 

Letting r = 0, shows $z)(-^) = I- Even if one assumes that ()19|) holds only when (J + r ■c?)/2 
is positive (that is when |r| < 1), then (fTU|) . combined with the linearity of completely 
determines ^d- Letting r = (1,0,0), we obtain $£i((/ + cri)/2) = (/ + (iiicri)/2. Hence, by 
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linearity, 



Similarly, $d(o"2) = '^220"2 and ^£,{(7-^) = d-^^a-^. Thus is a diagonal operator on A^2 with 
respect to the basis (ctq, cti, cr2, as). Of course, $d is completely determined by its action on 
this basis. 
Because 

crjCTjCTj = ±aj, 
the superoperator "if on 7^2, defined by 

3 

'^{p) = J2f3^aipai, (20) 

for some real constants will have J, cxi, ct2, and as eigenvectors. Thus \E' will 

equal if we can arrange to have \1/ yield the appropriate corresponding eigenvalues: 1, 
(ill, '^22; and (i33. This is a simple matter of solving the following linear system, the j-th 
equation of which is obtained by substituting p = aj into (j^Uj) : 

1 = /?o + /5i + /52 + /?3 (21) 
c^ii = + Pi - f32 - 

d33 = Po - (3i - (32 + (3^ . 

Let C denote the 4x4 matrix of coefficients of the preceding system and observe that C 
is a symmetric matrix such that = 4/. This observation permits quick solution of the 
system: 

_ 1 + C?ii + C?22 + (^33 o _ 1 + dll - d22 - ^33 fr,r,\ 

Po ^ ,Pi- ^ , {^^) 

1 — dll + d22 — dss n 1 ~ dll — d22 + dss 

(32 - ^ ,P3 - ^ • 

Hence we see that the superoperator has operator-sum decomposition given by 

3 

^d{p) = ^(SiffipCTi, 
i=0 

with the constants (3j given by ()22j). By Corollarv 111.21 is completely positive if and 
only if /5j > for j = 0, 1,2, and 3. Thus we have arrived at our desired characterization 
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of diagonal quantum superoperators on M.2- Observe that the nonnegativity of the /3j's is 
equivalent to the Algoet-Fujiwara conditions Q- 

Combining our work on diagonal superoperators with the factorization (jl8p . we find an 
operator-sum decomposition of the unital superoperator $jv/ defined by 

3 

^Af(p) = 5^ A(f/B^^,t^A)p(f/Ba^f/A)^ (23) 
i=0 

where M = BDA is the factorization of Proposition IIII. 1} where Ua and Ub are the unitary 
matrices given by (fT7|) . and where the scalars (3i are defined by ()22|) . As discussed above 
the superoperator will be completely positive if and only if the scalars (3i leading each 
summand are nonnegative. In P, Theorem 1(1)], Landau and Streater show that every uni- 
tal quantum superoperator on M.2 is a convex combination of unitary maps. Observe that 
fl2H|l recaptures the Landau- Streater result, and says a bit more: every unital superoperator 
$ on M.2 that preserves both Hermiticity and trace is a linear combination of unitary su- 
peroperators: $(p) = J2'i=o f^i^iP^i where Yl^=oPi — ^ /^^ ^^'^^j ^'"^ '^^^ necessarily positive, 
scalars Pi. 



V. DIAGONAL QUANTUM OPERATORS ON THE BLOCK REGION 

Let S = {0,1,2,3} be the index set for the Pauli matrices (including ctq = /) and let 
Sq = 5" \ {(0, 0, . . . , 0)} be the Cartesian product of n copies of S with the zero n-tuple 
removed. We represent the state p of a ra-qubit system in Bloch form 

^(^+V2"-n2"-l)'E r.A.j, 

where {Aj}^^"^~^ consists of all (appropriately normalized) n-factor tensor products of the 
Pauli matrices, excluding / = uo ® ■ ■ ■ ® o"o: 

{>^i}T=i'^ = { ^j^^ ^h ®(yh®---® ■ (ii'is, . . . , in) G ^0 }• 

Observe that {Aj}^^'^~^ together with Aq := // constitutes an orthogonal basis for A^2'' 
such that {Xi\Xj) = 26ij, where (-I-) is the Hilbert-Schmidt inner product: {A\B) = tr(A"''i?). 
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A basis should be ordered and we will use the "dictionary" ordering: 

Ao = — ctq ® (70 ® ■ ■ ■ (S> (To ® (To, Ai = — ao (g> do ® ■ ■ ■ ® o"o (S> (Xi, 

A2 = , q-o ® (To ® ■ ■ • ® CTq (g) 0-2, A3 = , (To ® (To ® ■ ■ ■ ® (To ® (T3, 

A4 = ^ (Tq ® ■ ■ ■ ® (Tq (S) (Ti (g) (To, . . . , A22n_i = ^ (T3 (T3 ® ■ ■ • (g) (T3 (g) 0-3. 

V2"-i V2^ 

Note that the dictionary ordering is equivalent to that produced by ordering according to 

the size of the index sequence iii2 ■ ■ - in associated with (g • ■ ■ (g (Tj„, where iii2 ■ ■ - in is 

interpreted as the base 4 representation of a number. 

Our goal is to characterize those (2^" — 1) x (2^" — 1) diagonal matrices D, with real 

entries djj, along the diagonal, such that $d : Ai2" defined by 

(/ + cf . a)) = ^ (/ + cDf. a) (24) 



is a quantum operator. Just as in the single-qubit setting, (j24|) together with the linearity 
of $£) yields 

= / and $d(Aj) = djjXj for j = 1, . . . , 2^" - 1. 

Because preserves Hermiticity, we know that it has an operator-sum decomposition; 
moreover, we know that has {Ajj^^To"^ as eigenvectors. The one-qubit situation, analyzed 
in the preceding section, suggests that has the form 

= f^-^^^P^^ (25) 

1=0 

for some real constants {pjfj^^ . 
Using (fTT|) . it is easy to check that 

A,A,A, = ±^A, for j G {0, 1, 2, . . . , 2^" - 1} (26) 

so that $D, given by ()25p. will have the right eigenvectors. We need only arrange for to 
have the correct eigenvalues (namely, 1, di i, . . . , (i22n_i^22"-i)- 

We need a way to keep track of the signs that appear on the right of ()26|) . Let C be the 
4x4 matrix of I's and — I's defined by 

(Ti(Tj(Tj = Cijaj 
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so that C is the matrix of coefficients of the system ^2\\ji of the preceding section. Note that 
C has I's along its main diagonal and C/2 is a symmetric, orthogonal matrix. Moreover the 
"one-qubit" Algoet-Fujiwara conditions — (3j > for j = 0, 1, 2, 3 where the /3j's are given 
by (j22j) — are equivalent to the requirement that the column vector 

1 

1 dii 
C 

.^33. 

have nonnegative components. 

Now let F be the 16 x 16 "sign" matrix corresponding to (f^H|) in the n = 2 qubit situation: 
XiXjXi = fijXj/2. It's not difficult to see that F = C ^ C; for example, to find the "upper- 
left" 4x4 block of F, one calculates 

SO that cqqC is the upper-left block of F, which is the appropriate submatrix in the Kronecker 
product. Thus when n = 2, the matrix of coefficients of the /5j's in the 16 x 16 system that 
results when Aq throught A15 are subsituted into is |C ® C, and thus, because the 
inverse of |C ® C is |C CSC, the necessary and sufficient conditions for to be completely 
positive is that 

have nonnegative components. Moreover, these components are precisely the /5j's in (I25|l 
that lead to an operator-sum decomposition of ^d- 
In complete generality, we have the following. 

Theorem V.l (Algoet-Fujiwara Conditions for n-qubits) The diagonal linear super- 
operator $D : A^2" M.2" defined by !\24l is completely positive if and only if the column 



1 

di,i 
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vector 



2n+i 



1111 
1 1-1-1 
1-1 1-1 
1-1-1 1 



1 



has nonnegative components. Moreover, if denotes j-th component of this column vector 
{forj = 1, . . . , 2^" ), then has operator-sum decomposition 

j=0 
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